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Abstract 

A novel approach to investigation of correlation effects in the electronic struc- 
ture of magnetic crystals which takes into account a frequency dependence of 
the self energy (so called "LDA++ approach") is developed. The fluctuation 
exchange approximation is generalized to the spin-polarized multi-band case 
and its local version is proposed. As an example, we calculate the electronic 
quasiparticle spectrum of ferromagnetic iron. It is shown that the Fermi liq- 
uid description of the bands near the Fermi level is reasonable, while the 
quasiparticle states beyond approximately 1 eV range are strongly damped, 
in agreement with photoemission data. The result of the spin-polarized ther- 
moemission experiment is explained satisfactory. The problem of satellite 
structure is discussed. 

I. INTRODUCTION 

The description of correlation effects in electronic structure and magnetism of iron-group 
metals is still far from the final picture and attracts continuous interest (see, e.g., |]~§]] and 
Refs therein). Despite many attempts, the situation is still unclear both theoretically and 
experimentally. For example, there is no agreement on the presence of 5 eV satellite in 
photoemission spectrum of iron |5]||, and on the existence of local spin splitting above 
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Curie temperature of nickel 0. The experimental data on the absence of spin-polarization 
in the thermoemission from cesiated iron || are still not understood completely |9j . From the 
theoretical point of view, different approaches such as the second-order perturbation theory 
T0|,|3|j, the three-body Faddeev approximation JTTJ] , and the moment expansion method fT2| 



were used. Unfortunately, the conditions of applicability of these schemes are not clear. 
Recently we investigated different approximations to the LDA-based correlated electronic 
structure of crystals with a local self-energy [[HI ( so called "LDA++" approach) and argued 



that for moderately strong correlations (the case of iron-group metals) one of the most 
efficient approaches would be the fluctuation-exchange (FLEX) approximation of Bickers 
and Scalapino [|14j]. Here we generalize LDA++ approach |13[ to the spin-polarized case and 



present some results for the quasiparticle spectrum of ferromagnetic iron. 



II. SPIN-POLARIZED MULTI-BAND FLEX APPROXIMATION 



Let us start with the general many-body Hamiltonian for crystal in LDA+U scheme ||15|| : 



H = H t + H v (1) 

AA'cr 

H u = \ ( A i A 2 M KK) ctacta'Cy^cy^, 

Z {\i}aa> 

where A = im are the site number (i) and orbital (m) quantum numbers; o =j, j, is the 
spin projection; c + ,c are the Fermi creation and annihilation operators; H t is the effective 
single particle Hamiltonian from LDA, corrected for double-counting of average interactions 
among correlated electrons |15|JT3|1 , and the Coulomb matrix elements are defined in the 
standard way 

(12 \v\ 34) = J drdrV*(r)V>2(r> (r - r) ^ 3 (r)^ 4 (r), (2) 



where we define for briefness Ai = 1 etc. Following Bickers and Scalapino |TJ| we introduce 
the pairwise operators corresponding to different channels, namely, particle-hole density: 



V2 

particle-hole magnetic channel: 



du - —7= (c ] t T c 2T + C^C 2 J. 



m 12 = C 1T C 2| 



m 12 = c u c 2T 

particle-particle singlet channel: 



S12 = -j= (c u c 2T - Ci T c 2i ) 



Sl2 — ( C 1T C 2| C l| C 2Tj > 

and particle-particle triplet channel: 

<ia = ( C U C 2T + c iT c 2|) 



*12 — ( C 1T C 21 + C U C 2T 

*12 = C 1T,I C 2J.,T 



l l2 — C 1T,I C 2J,,T • 



The bare vertex matrices corresponding to the different channels are defined as: 

Uf 2M = 2 (13 M 24) - (13|u|42) (3) 

^,34 = - (13 H 42) 

^12,34 = |[<12|^| 34) + (12 | V | 43)] 

£4,34 = \ [(12 34)- (12 | V | 43)]. 

Now we can rewrite the interaction part of Hamiltonian in Eq. ([!]) in two equivalent forms, 
corresponding to the particle-hole and the particle-particle channels: 



H v = -Tr ld*U d *d+ J2 m a * U m * m~ 

- ' Q = 0,± 



H v = -Tr I s * U s * s + t a *U t *t° 

- I a=0,± 



where (*) means matrix product, e.g. 



U 34 

( d * ud ) 12 =E^ d 



34,12 



Then we may repeat the usual derivation of FLEX equations for single band [1^] and 
multi-band [ P^| , P^[ cases, taking into account the spin-dependence of the Green function 
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T T c la (r) c+ (0) 



For the finite temperature (T> 0) FLEX equations has a "local form" in the Matsubara 
frequency (iu n ) or imaginary time (r) space (wherwheree u n = (2n + l)vrT, n = 0, ±1, ...) 
and it is very efficient to use the fast-Fourier transforms (FFT) with periodic boundary 
condition [|TJ],[nj. Time-frequency spaces are connected by 



e^Gl x ,(r)dr 



XX' 



T 



G a xx ,(iuj n ) 



and we will try to keep this dual (iu) — (r) notation to stress the numerical implementation 
of this LDA++ scheme. 

An unusual feature of the spin-polarized multi-band FLEX scheme is mixing of m° and 
(^-channels (as well as s and t°-channels) which related with non-zero value of such correlator, 
for example: 

((di2rnh)) c = \ (((c+c 2T + c+c 2i ) (c+c 4T - c+c 4i ))) c = -\ {gIgI - G l 23 G{x} , (4) 

where "c" denotes a connected part of the correlator (cf. |H[). As a result, we have the 
following expression for effective transverse susceptibility matrix: 



~i + u m *r tl \ 



lUJ) 



lUJ) 



(5) 



where 



(r) = -G° 23 (r) G£ (-r) 



(6) 



is an "empty loop" and F(iu) is its Fourier transform. The corresponding longitudinal 
susceptibility matrix has a more complicated form: 



where we introduce the supermatrix of p-h interactions: 



(7) 



U, 



( U d ' 



ph 



V 



u m 

and the matrix of bare longitudinal susceptibility: 

/ 



1 



Xo 



pTT _|_ pU pTT _ pU 



V 



(8) 



pTT _ pU pTT _|_ pU 

in the dd-, dm -, m°d-, and m°m°- channels (d,m° = 1,2 in the supermatrix indices). 
Similarly for the p-p channel we have: 



(9) 



where the supermatrix of p-p interactions is defined as: 



U. 



vv 



( U s ^ 



U* 



and the bare p-p susceptibilities are 

1 



12,34 2 



G\i (r) G23 (r) - G\i (r) (r) 



(10) 
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In the FLEX approximation we can calculate the electronic self-energy in terms of effec- 
tive interactions in various channels: 



where the Hartree-Fock contribution is equal to: 

£f£ = E (13M24)£<-(13M42K 4 

34 L a' 



(12) 



(13) 



with the occupation matrix n1 2 = G2i( T ~* ~ 0)i this contribution to X is equivalent to 
spin-polarized "rotationally-invariant" LDA+U method JHJ . 

The second-order contribution in the spin-polarized case reads: 



J; (13 H 74) G^(r) 

{3-8} 



85 | V | 26) £G&(t)G£(-<t) - (85 | v | 62) ^ 3 (r)^ 5 (-r) 
and the higher-order particle-hole contribution 



Sffi (r) 



E ^3% to ^34 (T) , 
34,ct' 



(15) 



with p-h fluctuation potential matrix: 

W aa ' (iu) = 



(iu) W n (iuj) 
W l Uiuj) WU(iw) 



(16) 



where the spin-dependent effective potentials are defined as 

1 



W 



n 



- 2 {u'* 
+u d * 



X M - X d o d 



*U d + U m * [ X mm - Xo m ] * u r 



dm dm 

A. A,0 



* u m + u m * 



md md 



u d } 
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W U = lm<i, 
-U d * 



X M ~ Xo 



* U d + U m * [ X mm - Xo m ] * u % 



dm dm 
X AO 



+ jjm _ jjm # 



md md 
X AO 



*U d } 



w u = u m * 



X - Xo 



*u r - 



= u m * 



X -Xo 



* u r 



Finally the higher-order particle-particle contribution corresponds to: 



with p-p fluctuation potential matrix: 



^42 (r) GJ, (r) , 



T aa> (iu) 



defined in terms of the spin-dependent p-p effective potentials: 



(17) 



T n = U S * [R ss - R S S ] *U S + UU 



pit T-)tt 

it — rtn 



*u r 



+u s * 



r>st r>st 
It — xt 



T)tS J-)tS 
It Xtg 



* u s 



T u = U S * [R ss - R S Q S ] *U S + U** 



It — Itn 



*u r 



-U s * 



T)St DSt 

it — JXq 



* W - W * 



jyts TDts 
It — JTq 



* u s , 



R + -R2 



* U l 



R — Rn 



* u*. 



Note that for both p-h and p-p channels the effective interactions, according to Eqs.( |l^ , |T7| ), 
are non-diagonal matrices in spin space, in contrast with any mean-field approximation like 
LSDA. 
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III. LOCAL APPROXIMATION 



The consideration of full non-local FLEX self-energy E^- (iu) (or momentum dependent 
Fourier transform E (k, iu)) in the framework of realistic multiband calculations is faced with 
significant computational difficulties ||16|| . On the other hand, the best local approximation 
which corresponds to the so called dynamical mean field theory (DMFT) |19| appeared to be 
very effective for the calculation of electron spectra of strongly correlated systems. Therefore 
it is reasonable to combine the FLEX approximation with DMFT rather than just neglecting 
momentum dependence in FLEX equations. In this respect we will take into account only 
on-site (Hubbard) interactions so all the bare vertices (Eq. ([!])) are diagonal in the site 
index i but are matrices in orbital indices {m}. We believe that this self-consistent local 
FLEX- approximation partly accounts for "vertex corrections" (due to the difference of G 
and Qi see bellow) which are absent in usual FLEX method. Moreover, it is known that 
even simple second-order approximation for self-energy combined with DMFT give rather 
good results for the description of correlated systems namely, Hubbard splitting, "Kondo 
resonances" etc. fl9f . 

The DMFT procedure on the lattice is as follows: one needs to find a self-consistent 
solution of functional equations 

Q- 1 (iu) = G- 1 (ioo) + S (iu;) (18) 

£ = £ [Go] , 

where 

G(iu) = J2[iu + A*-*(k) (19) 

k 

is the matrix (in orbital and spin indices) of the local Green function and £ (Eq.(|l^)) is 
our spin-polarized multi-band FLEX solution which depends on the effective media Green 
function Q Q . 

We could further reduce the computational procedure by neglecting dynamical interac- 
tion in the p-p channel since the most important fluctuations in itinerant electron magnets 
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are spin-fluctuations in the p-h channel. We take into account static ( of T- matrix type) 
renormalization of effective interactions replacing the bare matrix C/1234 = (12 \v\ 34) in Eqs. 



(14], Hip with the corresponding scattering T— matrix averaged over spins 



U = [l + U*U (iu = 0)]" 1 * U 



where Hn(iuj) is the Fourier transform of 



n? 2i34 (r) 



In the case of the single-band Hubbard model this approximation was found to be very 



reliable |[L8|| . The effects of the spin- dependence of the effective interaction matrix U will 



be considered elsewhere [ 2C . 



IV. COMPUTATIONAL RESULTS 

We have started from the spin-polarized LSDA band structure of ferromagnetic iron 
within the TB-LMTO method ]2l|] in the minimal s,p,d basis set and used numerical or- 
thogonalization to find the H t part of our starting Hamiltonian. We take into account 
Coulomb interactions only between d-states. The correct parameterization of the Hu part 
is indeed a serious problem. For example, first-principle estimations of average Coulomb in- 
teractions (U) lead to unreasonably large value of order of 5-6 eV in comparison with 
experimental values of the U-parameter in the range of 1-2 eV for iron ||. Semiempirical 



analysis of the appropriate interaction value [23| gives U ~ 2.3 eV. The difficulties with 



choosing the correct value of U are connected with complicated screening problems, defini- 
tions of orthogonal orbitals in the crystal, and contributions of the intersite interactions. In 
the quasiatomic (spherical) approximation the full [/-matrix for the d— shell is determined 
by the three parameters U, J and 5 J ||13| . Note that the value of intra-atomic (Hund) ex- 
change interaction J is not sensitive to the screening and approximately equals to 0.9 eV 
in different estimations P^] . We use the simplest way of estimating 5 J keeping the ratio 
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5 J/ J equal to its atomic value [p4|| . For the most important parameter U, which defines the 
bare vertex matrix (Eq. (|3])), we use the value U = 2.3 eV for most of our calculations and 
discuss the dependences of the density of states (DOS) as functions of U. To calculate the 
spectral functions 

A a (k,E) = --Tr L G a (k,E + iO) 

7T 

and DOS as their sum over the Brillouin zone we first made analytical continuation for the 
matrix self-energy from Matsubara frequencies to the real axis using the Pade approximation 
25 1 , and then numerically inverted the Green-function matrix as in Eq. (|19|) for each k-point. 



In the self-consistent solution of the FLEX equations we used 1024 Matsubara frequencies 
and the FFT- scheme with the energy cut-off at 100 eV. The sum over irreducible Brillouin 
zone have been made with 72 k-points for SCF-iterations and with 1661 k-points for the 
final total density of states. 

First we analyze the [/-dependence of electronic structure (Fig.l). Keeping J, 5 J fixed 
as described above we vary the average U parameter in the range from 2 to 6 eV. For com- 
putational simplicity a relatively high temperature value T = 1500 K was used. Note that 
the temperature in this approach is defined in terms of Matsubara frequencies. In contrast 



to the standard LDA calculations at the finite temperatures |26| we take into account the 
temperature dependence of the Bose degrees of freedoms through the p-h susceptibilities. It 
is known |27| that this source of temperature dependence is the most important for itinerant 



electron magnets. Although T is above T^ p (where T^ p = 1043 K is the experimental value 
of the Curie temperature of iron) all DOS (iVj (E) ,N± (E)) curves in the Fig.l show the 
spin splitting. In principle, one could calculate Tc using the temperature dependence of the 
uniform spin susceptibility ]TJ| but we have not done it yet. Nevertheless, we believe that 



this splitting in DMFT is a manifestation of the existence of local magnetic moment above 
Tq and is not connected with long-range magnetic order. This spin splitting is a charac- 
teristic feature of local approximations and reflects the Hubbard splitting in the one-band 



model [18[]. The local spin splitting above Curie temperature in iron is definitely observed in 
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many experiments, e.g. optics, photoemission etc. (see, e.g. |27j| and references therein). 



First of all, we see that the value of the total magnetic moment is weakly dependent on 
U and is of order of 2/^b for this temperature. The positions of the main peaks with respect 
to the Fermi energy roughly coincides with those in LSDA up to U = 2 eV. Starting from 
U = 3 eV satellites at E ~ —5 eV appears as well as additional many-body structure at 
E ~ 4 eV. Note that a weak satellite-like feature at E ~ —5 eV was observed experimentally 
in || although it was not found in ||. For U ~ 1 eV which is considered to be "experimental" 
value for iron H there are no any noticeable manifestations of this satellite. The boundary 
values of U of order of 2 eV with weak shoulder probably corresponds to the experimental 
situation in the best way. For U ~ 5 — 6 eV an essential part of the spectral density related 
to the many-particle peaks corresponding to the upper and lower Hubbard bands, which 
is unrealistic for such moderately correlated substances as iron. At U ~ 6 eV the empty 
quasiparticle minority-spin peak goes below the Fermi level which decreases the magnetic 
moment. 

The depolarization of states near the Fermi level is another important correlation effect. 
The decrease of the ratio P = [N r (E F ) - N l {E F )\ / [N r (E F ) + N x (E F )} is a typical sign of 
spin-polaron effects |{j[Z8|. In our approach this effects are taken into account through the 



W+l terms in the effective spin-polarized LDA++ potential. 

The energy dependence of self-energy in Fig. 2 shows characteristic features of moder- 
ately correlated systems. At low energies \E\ < 1 eV we see a typical Fermi-liquid behavior 
JmS (E) ~ — E 2 , dReTj (E) jdE < 0. At the same time, for the states beyond this interval 
within the ci-bands the damping is rather large (of the order of 1 eV) so these states corre- 
sponds to ill-defined quasiparticles, especially for occupied states. This is probably one of 
the most important conclusions of our calculations. Qualitatively it was already pointed out 



in Ref. ]TD[ on the basis of a model second-order perturbation theory calculations. We have 
shown that this is the case of realistic quasiparticle structure of iron with the reasonable 
value of Coulomb interaction parameter. 

Due to noticeable broadening of quasiparticle states the description of the computational 
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results in terms of effective band structure (determined, for example, from the maximum 
of spectral density) would be incomplete. We present on the Fig. 3 the full spectral density 
A a (k, E) including both coherent and incoherent parts as a function of k and E. We see that 
in general the maxima of the spectral density (dark regions) coincide with the experimentally 
obtained band structure. However, for occupied majority spin states at about -3 eV the 
distribution of the spectral density is rather broad and the description of this states in terms 
of the quasiparticle dispersion is problematic. This conclusion is in complete quantitative 
agreement with raw experimental data on angle-resolved spin-polarized photoemission |29| 
with the broad non-dispersive second peak in the spin-up spectral function around -3 eV. 



V. APPLICATIONS TO SPIN-POLARIZED THERMOEMISSION 

One of the most unexpected results concerning electronic structure of iron was obtained 
by spin-polarized thermoemission for cesiated iron ||. In this case the thermal current is 
determined by the states with the energy W = 1.37 eV above Fermi level which are in the 
region of quasiparticle DOS peak for minority spin. One could expect a strong negative 
spin polarization of the current (polarization ratio P as estimated from LSDA DOS is about 
-85% ). More accurate estimation which takes into account group velocities || results in P= 
-34% for the polarization ratio. Experimentally it found to be zero within the experimental 
error. 

To clarify the situation, we considered this effect on the basis of our LDA++ calculations. 
The current through the surface x = in the spectral representation is 

, ,dt(k) 

where dt (k) /9k = Vk is the group velocity operator (the matrix in orbital indices), f(E) 
is the Fermi distribution function. Taking into account only electrons moving towards the 
surface (k x > 0) and with the energy (E) above the barrier (W) the thermoemission current 
(cf. Ref || for the case of noninteracting electrons) could be expressed as: 



j x = Tr L* 



dt{k) 



dEf(E)A(k, E) 



(20) 
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oo , 

3l= E [dEf(E)Tr L A 



dt(k) 

dk T 



A(k,E) 



(21) 



k(k x >0) w 

Taking into account that T «W and averaging over the surface orientations (experiments 
H were carried out for polycrystalline samples), we found the following approximate formula 
for the polarization 



P 



(22) 



E T,Tr L 

i=x,y,z k 



dt(k) 



A a (k, W) 



This was calculated using numerical differentiation of t (k) matrix and summing up over 
1661 k points in the irreducible part of the Brillouin zone. We obtained P =-12% which is 
to compare with the value of —34% from LSDA calculations |§. The decrease of P is not 
a pure effect of damping of the quasiparticle states, but is the result of rather complicated 
cancellations of s,p,d -electron contributions. Therefore, one may conclude that there is 
no drastic discrepancy between experimental results || and theoretical description of the 
electronic structure of iron, inspite of the approximate character of our treatment of the 
thermoemission problem. For more accurate description one needs to consider , for exam- 
ple, an influence of the cesium layer and surface effects on the electronic structure of iron 
according to the experimental conditions. 



VI. CONCLUSIONS 

We have proposed a general scheme for investigation of the correlation effects in the 
quasiparticle band structure calculations for itinerant-electron magnets. This approach is 
based on the combination of the dynamical mean-field theory and the fluctuating exchange 
approximation. Application of LDA++ method gives an adequate description of the quasi- 
particle electronic structure for ferromagnetic iron. The main correlation effects in the 
electron energy spectrum are strong damping of the occupied states below 1 eV from the 
Fermi level Ep and essential depolarization of the states in the vicinity of Ep. We obtained 
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a reasonable agreement with different experimental spectral data (spin-polarized photo- and 
thermoemission) . The method is rather universal and can be applied for other magnetic 
systems, both ferro- and antiferromagnets. 
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FIGURES 




Energy, eV 

FIG. 1. Density of d-states and magnetic moments for ferromagnetic iron in the LSDA andthe 
LDA++ calculations for different average Coulomb interactions with J=0.9 eV and temperature 
T=1500 K. 



17 



-8 -4 4 8 

Energy, eV 



FIG. 2. Total density of states and d-part of self-energy for ferromagnetic iron with U=2.3 
eV and J=0.9 eV for the temperature T=750 K. Two different self-energies for t<2 g and e g d-states 
in the cubic crystal field symmetry are presented and four different lines corresponds to imaginary 
part spin-up (full line) and spin-down (dashed line) as well as real part spin-up (dashed-dot line) 
and spin-down (dashed-double-dot line). 
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FIG. 3. Spectral function of ferromagnetic iron for spin-up (a) and spin-down (b) and the two 
k-directions in the Brillouin zone compare with the experimental angle resolved photoemission and 
de Haas - van Alphen (at the Ep=0) points (from Ref. 3 ). 
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